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Abstract 



Using a recently proposed solution for an open antiferromagnetic spin-1/2 XXZ 
quantum spin chain with N (even) spins and two arbitrary boundary parameters at 
roots of unity, we compute the boundary scattering amplitudes for one-hole states. We 
also deduce the relations between the lattice boundary parameters appearing in the 
spin-chain Hamiltonian and the IR (infrared) parameters that appear in the boundary 
sine-Gordon S matrix. 



1 Introduction 



Factorizable 5* matrix is an important object of integrable field tfieories and integrable quan- 
tum spin cliains. As for the "bulk" case where the S matrix is determined in terms of 
two-particle scattering amplitudes, the "boundary" case can equally well be formulated in 
terms of an analogous "one-particle boundary-reflection" amplitude. These bulk and bound- 
ary amplitudes are required to satisfy Yang-Baxter [I]-[3] and the boundary Yang-Baxter 
[H [5] equations respectively. Methods based on Bethe equations have long been used to com- 
pute bulk two-particle S matrices [6]- [8]. In [8], Fadeev and Takhtajan studied scattering of 
spinous for the periodic XXX chain for both the ferromagnetic and antiferromagnetic cases. 
The bulk two-particle S matrix for the latter case coincides with the bulk S matrix for the 
sine-Gordon model [2] in the limit jS"^ —>■ Svr, where (3 is the sine-Gordon coupling constant. 
Much work has also been done on the subject for open spin chains [9]-[16] as well as for 
integrable field theories with boundary [5l|9]. In |i5j, Ghoshal and Zamolodchikov presented 
a precise formulation of the concept of boundary S matrix for 1 + 1 dimensional quantum 
field theory with boundaries such as Ising field theory with boundary magnetic field and 
boundary sine-Gordon model. For the latter model, the authors used a bootstrap approach 
to compute the boundary S matrix. They determined the scalar factor up to a CDD-type of 
ambiguity. Nonhnear integral equation (NLIE) [171 IE] approach has also been used to study 
excitations in integrable quantum field theories such as the sine-Gordon model [19]-[22] and 
open quantum spin- 1/2 XXZ spin chains [I2]^[15]. In fact, in [15], NLIE approach is used to 
compute boundary 5" matrix for the open spin-1/2 XXZ spin chain with nondiagonal bound- 
ary terms, where the boundary parameters obey certain constraint. The bulk anisoptopy 
parameter however is taken to be arbitrary. 

In this paper, we compute the eigenvalues of the boundary S matrix for a special case of 
an open spin-1/2 XXZ spin chain with nondiagonal boundary terms with two independent 
boundary parameters (with no constraint) at roots of unity, using the solution obtained 
recently [231 121] • The motivation for the performed computation is the fact that the Bethe 
Ansatz equation for this model is unchanged under sign reversal of the boundary parameters. 
Hence, the usual trick of obtaining the second eigenvalue of the boundary S matrix of an 
open spin-1/2 XXZ spin chain by exploiting the change in Bethe Ansatz equation under such 
sign reversal of the boundary parameters [lOl [TTl [151 IIS] would not work here. Consequently, 
identifications of separate one-hole states are necessary here. As far as the formalism goes, we 
follow the approach used earlier for diagonal open spin chains [TOlfTT] . This is a generalization 
of the method developed by Korepin, Andrei and Destri [HI [7] for computing bulk S matrix. 
The quantization condition discussed by Fendley and Saleur [9] is a crucial step for the 
calculation. The solution utilized here was derived for certain values of bulk anisotropy 



1 



parameter, /i in the repulsive regime (yU = G (0, |]) for odd p values. Hence, we focus 
only on the critical and repulsive regime, which corresponds in the sine-Gordon model to 
e [Atc, Svr) 0. One-hole excitations for this model occur in even N sector |25j in contrast 
to the diagonal open spin-1/2 XXZ spin chain where such excitations appear in the odd N 
sector [TTI. 



The outline of the paper is as follows. In Section 2, we briefly review the model. Previ- 
ously found Bethe Ansatz solution for the model is presented here [231 [21]. We also review 
the string hypothesis for two one-hole states. In Section 3, we proceed with the compu- 
tation of the scattering amplitudes. Since the Bethe roots for the model consist of "sea" 
roots and "extra" roots, we rely on a conjectured relation between the "extra" roots and 
the hole rapidity, which is confirmed numerically for system up to about 60 sites. We find 
that the eigenvalue derived for the open XXZ spin chain agrees with one of the eigenvalues 
of Ghoshal-Zamolodchikov's boundary S matrix for the one boundary sine-Gordon model, 
provided the lattice boundary parameters that appear in the spin chain Hamiltonian and 
the IR parameters that appear in Ghoshal-Zamolodchikov's boundary iS" matrix [5] obey the 
same relation as in [15] Q The problem of finding the second eigenvalue of the boundary S 
matrix requires the identification of an independent one-hole state. In contrast to previous 



studies [TOl 
parameters 



nl [T5l [T6] . where such state was found by reversing the signs of the boundary 
% similar strategy does not work here. Reversing the signs of the boundary 
parameters in the present case leaves the Bethe equation unchanged, hence giving the same 
one-hole state. Interestingly, a separate one- hole state with 2-string is found p5]. Using 
a conjectured relation between "extra" roots, hole rapidity and the boundary parameters, 
which is again confirmed numerically for system up to about 60 sites, we derive the remain- 
ing eigenvalue which also agrees with Ghoshal-Zamolodchikov's result. Finally, we conclude 
the paper with a brief discussion and possible future work on the subject in Section 4. 



2 Bethe Ansatz and string hypothesis 

We begin this section by reviewing recently proposed Bethe Ansatz solution [23l [24] for the 
following model [51 [2S] 

H = Tio +- sinh?7( cosech Q;_(Ti^ + cosech a_|_cr^) , (2.1) 



V'=8(7r-M) 

^Very recently, similar relations were found for the open XXZ spin chain with diagonal-nondiagonal 
boundary terms in [16j . 

•^In fact, there is a change £,± — > — f± in the Bethe equation for the diagonal case [TTj . 
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where the "bulk" Hamiltonian is given by 



Af-l 



^0 = ^ 5Z + <^n+l + cosh T] . (2.2) 



n=l 



In the above expressions, a^, a^, are the usual Pauli matrices, r/ is the bulk anisotropy 
parameter (taking values 77 = with p odd), a± are the boundary parameters, and is 
the number of spins/sites. Note that this model has only two boundary parameters. The 
most general integrable boundary terms contain six boundary parameters. In the present 
case, four other boundary parameters have been set to zero. We restrict the values of the 
remaining parameters, a± to be pure imaginary to ensure the hermiticity of the Hamiltonian 
(12. ip . The Bethe Ansatz equations are given by 



5{uf-r]) h^^){uf ) ~ Q2{uf+Tl) 



J = 1 , 2 , . . . , Ml , 

j = l,2,... ,M2. (2.3) 



where 



4/1 \\2Ar sinh2u sinh(2M + 4?]) 

o(u) = 2 (smhusmh(u + 2r7)) — — r — smh(u + + a_) 

^ ^ ^ V sinh(2M + 77) sinh(2M + 37]) ^ ' ^ 

sinh(u + 7] — a_) sinh(M + rj + q;_|_) sinh(M + rj — q;_|_) cosh^(M + rj) , (2.4) 

/.«(«) = 8 sinh^^-^^(n + 2,) cosh^(^ + ,) cosh(^ + 2^) ;,(2)(^^ ^ _ 2^p_5) 

sinh(2ri + Sr]) 

and 

Ma 

Qa{u) = Ylsmh{u-uf'^)siBh{u + uf'^ +ri) , a = 1,2, (2.6) 
i=i 

Ml and M2 are the number of Bethe roots, u^^^ and u^^^ (zeros of Qi{u) and Q2{u) respec- 
tively) . 



2.1 One- hole state 

In order to compute the spinon boundary scattering amplitude, we consider a one-hole state. 
The roots distribution for such a state was found in [2^. One-hole excitations for the open 
XXZ spin chain we study here appear in the even N sector. Hence, it is sufficient to review 
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following form 

Mf'^ + f, : J = l,2,... ,M(,,2) 



the results for even N case. The shifted Bethe roots Uj = Uj + ^ for this state have the 



f^K : j = l,2,...,M(ai) , , 

ia,2) . J ' a = l,2, (2.7) 



where /i = ^ and aJ-"'^^ are real. Here, M(i,i) = M(2,i) = f , and M(i,2) = 2±i, M(2,2) = 
The nX^""'^^ are the zeros of Qa{u) that form real sea ("sea" roots) and /iA^f'^'' are real parts 
of the "extra" roots (also zeros of Qa{u)) which are not part of the "seas". Hence, there are 
two "seas" of real roots. We employ notations used in 



sinh (/i(A - f )) ^ 2fi' cosh (/i(A - f )) ^ ^ 

Rewriting bulk and boundary parameters [13], rj = ifi, a± = i/iaj^ the Bethe Ansatz 
equations (12. 3p for the "sea" roots then take the following form 

e,{Xf''^r'''-'[gi{Xf''^)e,^2aA^^^^^^ (2.9) 

N/2 (p-l)/2 

= - n - ").2(a;"' h- ai^^")] n - a?'=')*(a5"' + Af =') 

fc=l fc=l 
and 

,^(^(2,l))2A.+l^^(^(2,l))-l ^2.10) 

N/2 (p+l)/2 

= - n [-^(^f - ^1"V2(a;"' + At">)] n [*(Af - AL"')*(Af H- At'^=>; 

fc=l fc=l 

respectively, where j = 1 , . . . , y. These equations can be re-expressed in terms of counting 
functions, h'-'^(A) as 

h«(Af')) = J„ 1 = 1,2 (2.11) 

where h^')(A) are given by 

h^'^(A) = ^{(2iV + l)gi(A) -ri(A) -gi+2a_(A) -gi_2a_(A) -gi+2a+(A) -gi-2a+(A) 

N/2 (p-l)/2 
A;=l fc=l 



"'The string hypothesis ()2.7p holds true only for suitable values of a±, namely < \a±\ < -^^ii , 
a+a_ > , where v = p + 1 



4 



and 



h^'\X) = ^\{2N + l)q,{X)-n{X) 



271 

N/2 (p+l)/2 
A:=l k=l 

In the above equations, qn{X) and r„(A) are odd functions defined by 

Qn{X) = TT + z Inen(A) = 2 tan~^ (cot(?7./i/2) tanh(/iA)) , 

r„(A) = ^lnf7„(A). (2.14) 

Further, {Ji, J2, ■ ■ ■ , Ja} is a set of increasing positive integers that parametrize the state 
@. For states with no holes, the integers take consecutive values. For one-hole state, there is 
a break in the sequence, represented by a missing integer. This missing integer J, fixes the 
value of the hole rapidity. A, according to 

h«(A) = h(2)(A) = J. (2.15) 

If the hole is located to the right of the largest "sea" root (A^j^'^"*), then J = [h(')(oo) — 

T 

hW(A^^''^)J. See for more details. For later use, we next define the densities of "sea" 
roots as 

where 1 = 1,2 

The functions fl2.14p have the following derivatives which prove to be essential to the analysis 
in following sections. 

Id /i sin(ra/i) 

O-niX) = 7^TT?n(A) - 



27r dX rr cosh(2/iA) — cos(n/i) 

, /,x ^ d. u sin(n/i) ,^^r.^ 

bnW = ^TT^n(A) = -^ j ■ (2.17) 

ZTT d\ Ti cosh(2/iA) + cos(n/ij 

2.2 One-hole state with 2-string 

In addition to the one-hole state considered in last section, there is another one-hole state. 
This state is the only remaining one-hole state, which also has a 2-string. In this section. 



^In principle, there are two such sets of integers, {j/^^} and {jf'^^ corresponding to the two counting 
functions, h(^)(A) and h'^^(A) respectively. But, in fact these two sets of integers are identical. Hence we 
choose to drop the superscript, / from Jj in p. lip . 
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we give some brief information on the state. The shifted Bethe roots u^""^ 
state have the following form 



u^f* + ^ for this 



/iAf ''^ + f , 



1,2, 
1,2, 



(a,l) 
(a,2) 



1,2, 



(2.18) 



(a) »7 



where Ag"'' , ^ 



^ and Aj'^'^'' are real. Here, M(i i) 

_ PllS 1 — f — ,, 

yuA^,"'^"* are real parts of the "extra" roots (also zeros of Qa{u)) which are not part of the 
"seas". For this state, we also have /iAq^^ the real parts of additional "extra" roots that 
form a 2-string. 

The counting functions for this state are given by 
h^'^A) = ;^{(2iV+ l)gi(A) -ri(A) -gi+2a_(A) -gi-2a_(A) -gi+2aJA) -gi„2a-,(A) 



M, 



(2,1) 



f - 1, and M(i,2) 



p-i 

2 ' 



M, 



(2,2) = As before, fiXj are the zeros of Qaiu) that form real sea ("sea" roots) and 



and 



-E[^2(A-An+g2(A + Af^^) 
fc=i 

-g3(A - A^')) - g3(A + A^'^) - q,{\ - \f) - q,{\ + 
h(2)(A) = ^{(2iV+l)gi(A)-ri(A) 



(p-3)/2 



(2.19) 



(p-l)/2 



aM) + .,(a + aM; 



-g3(A - A«) - g3(A + A«) - q,{\ - A«) - gi(A + A^^^) 
The Bethe Ansatz equations for this state take the following form 

h«(Af^)) = J„ / = 1,2 



(2.20) 



(2.21) 



where {Ji, J2, . . . , Jk_i\ is a set of increasing positive integers that parametrize the state. 
The hole for this state breaks the sequence, represented by a missing integer. As before, the 
missing integer J, enables one to calculate the hole rapidity, A using 

h«(A) = h(2)(A) = J. (2.22) 

If the hole appears to the right of the largest "sea" root (A^'^"* ), then J = [h('^(oo) — 
h*^'^(A^'^''^)J . More on this state can be found in 



3 Boundary S matrix 



In this Section, we give the derivation for the boundary scattering amphtudes for one-hole 
states reviewed in Section [2j 



3.1 Eigenvalue for the one- hole state without 2-string 



First, we consider the state reviewed in Section [2711 From fl2.12p . (12.131) . (12.161) and (I2.17p . 
one can solve for the sum of the two densities. We recall the results below [21], 



PtotaliX) = P^'\X)+P^''>{X) 

= 4s(A) + li?+(A) 



(3.1) 



where s(X) — — r;/-^ 
by 



and -R+(A) is the inverse Fourier transform of R+iu}) § which is given 



1 + 02 (cj) 



2ai{uj) + 2a2(t^) - 2hi{uj) - ai+2a„(t^) - ai_2a_(t^) - ai+2a+(t^) - ai_2a+(c^) 



2 



P+1 
2 



262(cc')( J]cos(Af^^cj) + ^cos(Ap'^^cu)) +4a2(cu) cos(Au;) 



(3.2) 



and 



fe=i 



1=1 



sgn(n 



sinh ((z/ — \n\)uj/2) 



sinh {vuj/2) 
sinh {nuj/2) 



sinh (i/u;/2) 



< |n| < 2z/, 
< 3?e n < z/ . 



(3.3) 
(3.4) 



are the Fourier transforms of f[2.17p . The presence of "extra" roots, A^"'^"* and the hole 
rapidity. A, are to be noted herqj. Henceforth, we shall denote A[."'^'* simply as A^""*. Morever, 
momentum of the excitation is given by 



p{X) = tan ^ (sinh(7rA)j — — 



vr 



(3.5) 



®Our conventions are 



e'^^ fix) dX., 



/(A) 



2n 



-iojX 



/(w) dui . 



''Energy carried by the hole is given by E{X) — ^ 2" ^ cosIiVa) ' 



Such an expression for spinon was derived 



ini 
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From (13. Sp . one gets s(A) = Consequently, using fl2.16p . one rewrites (13. ip as 

1 dhtotai{X) 2rfp(A) 1 

where hto^^^A) = h(i)(A) + h(2)(A) and ^^^^^^ = Pto^KA)- After integrating ^ with 
respect to A, taking hmits of integration from to A, one finds □ 

Wotai{\) = h«(A) + h(2)(A) = -iVp(A)+ /" d\R+{\) (3.7) 

Since h'^^)(A) = h'^^)(A) G positive integer and -R+(A) is an even function of A, multiplying 
the resulting expression by 2m and exponentiating gives 

g2ip(A)iVgf /_\dAK+(A) ^ ^ 

Next, let us compare equation (13.81) to the Yang's quantization condition for a particle on 
an interval of length A^, 

g2ip(A)7V^^^. a_) I A, (±)) =1 A, (±)) (3.9) 

where -R(A; a±) are the non-diagonal boundary S matrices and | A, (±)) denote the two 
possible one-hole states. Note that the ± in | A, (±)) represents two posssible one-hole states 
and not the right and left boundaries. The expression e~ ^-'^''''^^^''^^ then, should be equal 
to one of the two eigenvalues of the Yang matrix Y{X) defined by 

Y{\) = R{X-a+)R{\-a^) (3.10) 

Defining this eigenvalue as q;(A, a+)a(A, a_), where -|- and — denote the right and left bound- 
aries respectively, (13. 8p can be rephrased as 

, a+)a;(A, a_) = 1 (3-11) 

The problem thus reduces to evaluating the following 

a(A, a+)a(A, a_) = e"^ ^-^ '^^^+(^) (3.12) 

After some manipulations, we have the following. 



oo 



duj 







a2{uj) 



a(A, a+)a(A, a_) = exp < 2 / — sinh(2iAci;) ^ H — — -[0,2(20;) + ai(2u;) — 6i(2u;) 



l + aaM l + a2(2uj) 



^(ai+2a_(2u;) + ai_2a_(2u;) + ai+2a+(2u;) + ai_2a+(2t^)) 



p— 1 P+1 
2 



S2(2cj)(J]cos(2Af^cj) + J]cos(2Ap^cu))] } (3.13) 



k=l 1=1 



^Since we are only able to determine the scattering amplitudes up to a rapidity-independent factor, the 
additive constant p(0) from the integration is ignored in p.7p . 
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Further, using fl3.3p and (13 ■4p . one gets 



a(A, a+)a;(A, a_) = exp <j 2 
+ 








■ sinh(2iAco') 



sinh(c(j) 



+ 



-2 sinh(3cj /2) sinh((z/ - 2)cj /2) 
. sinh(2cj)sinh((z/- l)cj/2) 

sinh((— z/ + 2a_ — l)c(j) 



sinh((z/ — 1)^^;) cosh(co') 2 sinh((z/ — l)a;) cosh(a;) 

sinh((z/ — 2a„ — l)ci;) 

+ ^T-TTT TV-\ t/t + ^^-^^+) 

2 sinn((z/ — l)uj) cosn(u;j 

sinh(c<j) 
sinh((z/ — l)^;) 



2 



+ 



^ cos(2Af + J] cos(2Ap)^))] } (3.14) 



where (a_ — a+) is a shorthand for two additional terms which are the same as the third and 
fourth terms in the integrand of (I3.14p . but with a_ replaced by a+. The integrals involving 
"extra" roots A^^"* and A^^^* yield 

exp (2 r — sinh(2au;) '^^^'^\ . (V cos(2A1')c.) + V cos(2A«a;)) 
Wo ^ smh((z/-lV)'^ 



p— 1 P+1 
2 2 



nnt//(^"'.^!".w(Af,A,'".A+^' 



(3.15) 



where /i' 



' - ^ and 

I/— 1 



/(Ar,A«,A) 



sinh(^(A + Al^) + |- 



.))cosh(^(A-Af + f + ^)) 



sinh (f (A + Al^) - f + ^)) cosh (f (A - A^ - f - 



t)) 



X 



sinh (^(A + Ap) + t + ^)) cosh (^(A - A^^^ + f - f.)) 



sinh (^(A + A« - t - ^)) cosh (f (A - Aj^^ - | + f.)) 
After evaluating the rest of the integrals, (I3.14p becomes 



(3.16) 



p— 1 P+1 
2 2 



a(A, a+)a(A, a_) = 5o(A)25i(A, a^)S^{X, «+) H H ^!'^ >^)fi>^k\ >^l'\ A + 



A:=l /=1 



(3.17) 



where 



■ COS. 



vr 



h(/.'A) n 



r[^(4n + l-2a)" 


r 


^(4n + 3-2zA) + l 


r[^(4n + l + 2a) 


r 


^(4n + 3 + 2a) + 1 



X 



X 



r 


^(4n + 4 + 2iA) 


r 




r 


^(4n + 4-2a) 


r 


^(4n-2zA) + 1 



r2 


^^(2n-^A) + i 


r2 


^(2n + l + zA) + i 




r2 


'^(2n + 2 + a) + l 


Y2 


_i_(2n + l-zA) + i_ 



(3.18) 



5'i(A,a±) 



- Jcosh(/i'(A + ^(z/ - 2a±))) cosh(/i'(A - ^(z/ - 2a±))) 

TT V 2 Z 



X 



X 



^ r [^(2n + 1 + a - - 2a±)) + ij r [^(2n + 1 + zA + i(z/ - 2a±)) + i 

n=0 



r I ^(2n + 1 - zA - i(z/ - 2a±)) + i T -l_(2n + 1 - zA + i(z/ - 2a±)) + i 

S.19) 



r 


_i_(2n-a-i(z/-2a±)) + i 


r 


_i_(2n-zA + i(z/-2a±)) + i 




r 


_l_(2n + 2 + ^A-i(z/-2a±)) + i 


r 


^{2n + 2 + zA + i(z/ - 2a±)) + | 



The values of the "extra" roots are dependent on the hole rapidity, A and the boundary- 
parameters, a±. Hence, it is sensible to expect a relation between these "extra" roots, 
jA^^-*, A^^-*!, the boundary parameters, a± and the hole rapidity, A. Consequently, one needs 

to express the right hand side of (13.151) in terms of purely a± and A to complete the derivation. 
To look for this additional relation, we begin with the information contained in the difference 
of the two densities, p^^-'(A) — p*^^)(A). This leads to the following. 



(3.20) 



where R-{X) has the following Fourier transform 

1 



Rjuj) 



ai+2a_('^) - ai-2a_(t^) " ai+2a+('^) " ai-2a+ (t^) 



1 - 0-2(^1;) 

- 262(.;)(X]cos(a1^)^) - 5^cos(Ap)^)) 



2 



k=l 



Analogous to (13. 6p ) one gets 



1 dhdiffjX) 
N dX 



N 



R-W 



(3.21) 



(3.22) 



where KiffiX) = h^^\X) - h(2)(A) and ^^^^^^ = Pdiff{X). Further, integrating (Km with 
respect to A, taking limits of integration from to A as before, one finds 



W(A) = h«(A)-h(2)(A)= / dAi?_(A) 

'0 



(3.23) 
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Since h(^)(A) = h(^)(A) G positive integer, using the fact that -R-(A) is an even function of A 
and exponentiating fl3.23p we get 



2 2 



^(A,a+)^(A,a_) JJ JJ 



/(Ai^A«,A) 



\i=i\ /(Al^\ Ap\ A + 



(3.24) 



where g{X, a±) = ^"^^,^1/'''^ 2a±))+^smh(/J^A) ^ Next, an important observation is the following 



cosh{ ^ {u-2a±))-i smh(/i' A) 



relation (as — > ool 



2 2 



nn/(^?^^j'^^)=-i (3-25) 

fc=i «=i 

for which we provide numerical support in Table 1. Although the results shown in Table 1 
are computed for the case where the hole appears to the right of the largest "sea" root, we 
find similar results for other hole locations. From (13.241) and (13.251) . it also follows that 



2 2 



k=l 1=1 



(3.26) 



N 


nEln3i/(Af,A«,A),p = 3 


nEn3i/(Af,A«,A),p = 5 


24 


-0.999364 + 0.0356655 i 


-0.999334 + 0.036496 i 


32 


-0.999421 + 0.0340133 i 


-0.999333 + 0.036522 i 


40 


-0.999466 + 0.0326686 i 


-0.999334 + 0.036486 i 


48 


-0.999502 + 0.0315413 i 


-0.999337 + 0.036419 i 


56 


-0.999532 + 0.0305749 i 


-0.999340 + 0.036336 i 


64 


-0.999558 + 0.0297318 i 


-0.999343 + 0.036243 i 



Table 1: Ilii IlSi /(A1^\aJ'\A) for p = 3 (a+ = 2.1, a_ = 1.6) and p = 5 (a+ = 3.3, 
a_ = 2.7), from numerical solutions based on = 24 ,32 ,. . . ,64. 

We stress here that the values of A^^'' and A^^^* used in computations above strictly satisfy 



the Bethe equations (12.31) . Finally, we can rewrite (13.171) as 

q;(A, a+)a(A, a_) = S'o(A)^S'i(A, a_)S'i(A, a+)(7(A, a+)5'(A, a_ 



(3.27) 



up to a rapidity-independent phase factor. Subsequently, the complete expression for each 
boundary's scattering amplitude is given by (up to a rapidity-independent phase factor) 



a(A,a±) = 5'o(A)5'i(A,a±)5f(A,a±) 
where -|- and — again denotes right and left boundaries respectively. 



(3.28) 
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3.2 Eigenvalue for the one-hole state with 2-string 



We now consider the one- hole state with a 2-string, reviewed in Section [521 The computation 
of the eigenvalue for this state is identical to the one given above. Hence, we skip the details 
and present the result. Analogous to (13.141) . we have 

-2sinh(3u;/2)sinh((i/ - 2)uj/2) 

to 

sinh(ti;) 



/3(A, a+)/5(A, a_ 



exp s 2 / — smh.[2iXiu) 



+ 



+ 



. sinh(2tu)sinh((z/- l)cj/2) 
sinh((— z/ + 2a_ — 1)uj) 



sinh((z/ — l)^^) cosh(ti;) 2 sinh((z/ — l)u;) cosh(co') 



^ sinh((z/ — 2a_ — l)c(j) 
2 sinh((z/ — l)u;) cosh(co') 

P-3 

sinh(u;) ^v^„„„/o\(2) 



£-1 

2 



+ 



sinh((z/ — 

which after evaluating the integrals yields 



sinh((z/- IV)^^ V z 

'^^^^^'^-L(cosh(2zA(^)^) + cosh(2^Ai^).;))l ) 



(3.29) 



P{X, a+)/3(A, a_ 



5o(A)^5i(A, a_)5i(A, a+)w{X^^\ X)w{X^'> , A) 



p—3 p— 1 
2 2 



X 



(3.30) 



fc=l z=i 



where 



w{Xl,^\X) 



cosh(/i^(A + Aj)"^ + i)) cosh(/iXA - A^°^ + z)) 
^ cosh(/i'(A + X^o^ - i)) cosh(/i'(A - X^o^ - i)) 



(a) 



a = l,2. 



(3.31) 



As before, (a_ — ^ a+) represents two additional terms which are the same as the third and 
fourth terms in the integrand of (13.291) . but with a_ replaced by a+. We proceed to make 
the following conjecture (as oo) to complete the derivation. 



p—3 p—1 
2 2 



W 



(A«, A)^(A?\ A) n n A)/(Ar\ A«, A + ^) = ^(A + ^, a,),(A + ^, a_) 

^ jjj jJi 



k=l 1=1 



(3.32) 

Like (13.251) , we provide numerical support for (I3.32p in Table |2] where we compute the ratio 
(f) = where di and ^2 are the left hand side and the right hand side of (13.321) respectively, 
for systems up to 64 sites. We believe this supports the validity of (13.321) at A^ — > cxd. The 
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values of A^^"* , A^^'' , Aq^'' and Ao^"" used in computations are obtained by solving numerically 
the Bethe equations (12.191) and (12.201) for the "sea" roots and (12.31) for the "extra" roots. The 
correctness and validity of such numerical solutions are checked by comparing them with the 
ones obtained from McCoy's method for smaller number of sites, e.g., N = 2 ,4 and 6^. We 
stress here that although the results obtained in Table [2] are computed for J = 1, namely 
the case where the hole appears close to the origin, similar results are found for other hole 
locations, e.g., J = 2 , 3 , . . .. 



N 


0, p = 3 


0, p = 5 


24 


0.967073 + 0.254500 i 


0.990295 + 0.138982 


i 


32 


0.981063 + 0.193688 i 


0.994434 + 0.105361 


i 


40 


0.987716 + 0.156259 i 


0.996308 + 0.085849 


i 


48 


0.991392 + 0.130928 i 


0.997674 + 0.068166 


i 


56 


0.993634 + 0.112654 i 


0.998065 + 0.062174 


i 


64 


0.995102 + 0.098852 i 


0.998407 + 0.056428 


i 



Table 2: for p = 3 (a+ = 2.1, a_ = 1.6) and p = 5 (a+ = 3.2, a_ = 2.7), from numerical 



solutions based on = 24 ,32 ,. . . ,64. 



Using (13.321) . the other eigenvalue for the Yang matrix (I3.10p becomes 
/?(A,a+)/3(A,a„) 



5'o(A)^S'i(A, a+)Si{\, a^)g{\ + — , a+)g{\ + — , a_ 



(3.33) 



hence giving the following for each boundary's scattering amplitude (up to a rapidity- 
independent phase factor). 



/3(A, a±) = 5'o(A)5'i(A, a±)g{X + — , a±) 



(3.34) 



3.3 Relation to boundary sine-Gordon model 

Next, we briefly review Ghoshal-Zamolodchikov's results for the one boundary sine-Gordon 
theory [5]. We borrow conventions used in flbl [27]. Ghoshal-Zamolodchikov's results imply 
that the right and left boundary S matrices R{9 ; ri±, j^) are given by 



R{e ; rj, ^, 7) = ro{9) n {6 ; r^, M{e ; 77, ^, 7) 
where M has matrix elements 

mil 



17121 rn22 



(3.35) 



(3.36) 



^We are only able to use McCoy's method to exactly solve for the Bethe roots for systems up to only 6 
sites due to computer limitations. 
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where (?7±, 'd±, 7±) are the Ghoshal-Zamolodchikov's IR parameters and 6 is the hole-rapidity. 
Further, 

mil = cos T] cosh i} cosh(r0) + i sin t] sinh sinh(r^) , 
777,22 = COS ?7 cosh i} cosh(r0) — i sin rj sinh sinh(r6') , 
mi2 = ie*''' sinh(r0) cosh(r0) , 

m2i = ie~^^ sinh(re) cosh(r^) . (3.37) 

where r = is the bulk coupling constant. The scalar factors have the following integral 
representations [15l [27] 

r dco ,^ , ,sinh((z/-2)cj/2)sinh(3cj/2) 

^\ Jo uj ^ ' ' sinh((z/- lV/2)sinh(2tu) 

ri(^;r/,^) = L_^^(^, ^) ^(,^, ^) , (3.38) 

cos ri cosh V 

where 



oo 



, , du ,, ^, ,, ,^ ,, cosh((i/ — l)u;x/7r) ,, 

a(.,9)=e.p^2 / - sm((,. - g)W(2.)) ..^B. ! (2.)) ^^^^^^ ^^>^^^^ p.39) 







Our result f l3.27p and (13.331) agree with the eigenvalues of -R(6' ; ?7+, -(9+, 7+)-R(6' ; ?7_, -i?-, 7_), 
provided we make the following identification, 

u' 

V± = y(i^-2a±) 
^ = ttA (3.40) 

In addition to (13.401) . one should also take '(9± = 7± = 0, since they are related to the 
lattice parameters that appear in the spin chain Hamiltonian, (12.11) which have been set to 
zero. Refer to the discussion following (12.21) . The same expression is given in [15j for the 
corresponding open XXZ spin chain with nondiagonal boundary terms but with a constraint 
among the boundary parameters, hence suggesting that (I3.40p holds true in general. As 
noted above, the eigenvalues, (13.281) and (I3.34p agree with the sine-Gordon boundary S 
matrix eigenvalues. Hence the two eigenvalues can be related as follows, 

a{X, a±) cosh(^(z/ — 2a±)) + i sinh(yu'A) 



P{X,a±) cosh(^(i/ - 2a±)) - i sinh(/i'A) 

4 Discussion 



(3.41) 



Based on a recently proposed Bethe ansatz solution for an open spin- 1/2 XXZ spin chain with 
nondiagonal boundary terms, we have derived the boundary scattering amplitude (equation 
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fl3.28p ) for a certain one- hole state. We used a conjectured relation between the "extra" 
roots and the hole rapidity, namely (13.251) . which we verified numerically. This result agrees 
with the corresponding S matrix result for the one boundary sine-Gordon model derived by 
Ghoshal and Zamolodchikov [5] , provided the lattice and IR parameters are related according 
to (13.401) . We obtained the second eigenvalue (13.341) by considering an independent one-hole 
state with a 2-string. This scattering amplitude, derived for the one-hole state with 2- 
string also agrees with Ghoshal-Zamolodchikov's result following conjecture (13.321) . which 
we verified numerically and identification (I3.40p . It would be interesting to derive (I3.25P and 
(I3.32P analytically. 

It will also be interesting to study the excitations for the more general case of the open 
XXZ spin chain, namely with six arbitrary boundary parameters and arbitrary anisotropy 
parameter, and derive its corresponding S matrix. Solutions (spectrums) have been pro- 
posed for the general case, using the representation theory of q-Onsager algebra [28] and 
the algebraic-functional method [29]. However, Bethe Ansatz solution for this general case 
has not been found so far although such a solution has been proposed lately for the XXZ 
spin chain with six boundary parameters at roots of unity [30]. In addition to the bulk 
excitations, one can equally well look at boundary excitations although this can be rather 
challenging even for the simpler case of spin chains with diagonal boundary terms. It is 
therefore our hope that some of these issues are addressed in future publications. 
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